LUMINOSITY AND STELLAR MASS FUNCTIONS OF 
LOCAL STAR-FORMING GALAXIES 



Pablo G. Perez- Gonzalez, Jesus Gallego, Jaime Zamorano 1 
Almudena Alonso-Herrero 2 
Armando Gil de Paz 3 
and 

Alfonso Aragon-Salamanca 4 
ABSTRACT 

We present the optical and near-infrared luminosity and mass functions of the 
local star-forming galaxies in the Universidad Complutense de Madrid (UCM) 
Survey. A bivariate method that explicitly deals with the Ha selection of the 
survey is used when estimating these functions. Total stellar masses have been 
calculated on a galaxy-by-galaxy basis taking into account differences in star 
formation histories. The main difference between the luminosity distributions 
of the UCM sample and the luminosity functions of the local galaxy population 
is a lower normalization ((f)*), indicating a lower global volume density of UCM 
galaxies. The typical near- infrared luminosity (L*) of local star- forming galaxies 
is fainter than that of normal galaxies. This is a direct consequence of the lower 
stellar masses of our objects. However, at optical wavelengths (B and r) the 
luminosity enhancement arising from the young stars leads to M* values that are 
similar to those of normal galaxies. The fraction of the total optical and near 
infrared luminosity density in the local Universe associated with star-forming 
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galaxies is 10-20%. Fitting the total stellar mass function using a Schechter 
parametrization we obtain a = —1.15 ± 0.15, log(.M*) = 10.82 ± 0.17 M & and 
log(0*) = — 3.04±0.20 Mpc~ 3 . This gives an integrated total stellar mass density 
of io 7 - 83±0 ' 07 .M© Mpc -3 in local star-forming galaxies (H = 70 kms -1 Mpc -1 , 
Qm = 0.3, A = 0.7). The volume-averaged burst strength of the UCM galaxies is 
b = 0.04 ± 0.01, defined as the ratio of the mass density of stars formed in recent 
bursts (age < 10 Myr) to the total stellar mass density in UCM galaxies. Finally, 
we derive that, in the local Universe, (13 ± 3)% of the total baryon mass density 
in the form of stars is associated with star-forming galaxies. 

Subject headings: galaxies: luminosity function, mass function — galaxies: star- 
burst — infrared: galaxies 

1. Introduction 

The star formation activity is a topic of growing interest as we extend the study of 
the galaxy population to higher and higher redshifts. Part of this interest is due to the 
fact that many distant objects show an intense starbursting activity (see, e.g., Ellis 1997; 
Ferguson et al. 2000, and references therein). One of the issues that is being addressed is the 
characterization of the global properties of the galaxy population as a function of lookback 
time. One of the most important of these properties is the comoving star formation rate 
(SFR) density of the Universe (see, for example, Haarsma et al. 2000; Somerville et al. 2001; 
Rowan- Robinson 2001; Lanzetta et al. 2002). Another issue of cosmological interest is the 
luminosity evolution of the galaxies along the Hubble time. One way to study this evolution 
is through the measurement of the number density of galaxies of a given luminosity as a 
function of redshift, i.e., the time-evolution of the luminosity function (LF). A considerable 
amount of effort has been devoted to the determination of the LF at a variety of redshifts 
and spectral ranges (among others, Efstathiou et al. 1988; Steidel et al. 1999; Loveday 2000; 
Kochanek et al. 2001; Cole et al. 2001; Norberg et al. 2002). 

Mass is even more fundamental than luminosity among the parameters governing the 
formation and evolution of galaxies. Indeed, it may be possible to distinguish between 
hierarchical and single collapse models of galaxy formation using the mass of the observed 
galactic structures and the history of their assembly as the differentiating parameters (see 
Aragon-Salamanca et al. 1998; Baugh et al. 1998; Fukugita et al. 1998; Kauffmann et al. 
1999). Therefore, it is extremely important to characterize the stellar and total masses of 
galaxies at all redshifts and the number density of objects of a given mass, i.e., the stellar 
mass function (SMF). 
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In this Letter, we have used the multi-wavelength dataset for the Universidad Com- 
plutense de Madrid (UCM) Survey of emission-line galaxies (ELGs, Zamorano et al. 1994, 
1996) to estimate the LFs and SMF of a complete sample of local star-forming galaxies in the 
optical and near-infrared (NIR). Given that the UCM Survey galaxies were selected by their 
active star-formation, our results are directly comparable to the ones achieved for samples 
of star-forming galaxies at high redshift, though one must be aware of the differences on the 
selection techniques. Unless otherwise indicated, throughout this Letter we use a cosmology 
with H = 70 kms" 1 Mpc" 1 , tt M = 0.3 and A = 0.7. 

2. The sample 

The UCM sample contains 191 local star-forming galaxies at redshifts lower than 0.045, 
with an average (z) = 0.026. These objects were selected from an objective-prism survey 
carried out with the Schmidt Telescope at Calar Alto Observatory (Almeria, Spain) and 
centered at the wavelength of the Ha nebular emission (Zamorano et al. 1994, 1996). Com- 
parisons with other samples of star-forming galaxies selected using a variety of techniques 
reveal that the UCM sample is representative of the local population of star-forming galaxies 
(Gallego et al. 1997; Salzer et al. 2001; Salzer 1989; Treyer et al. 1998). Indeed, the SFR 
density of the local Universe derived from the UCM sample is reasonably consistent with 
that derived from other samples (see, e.g., Treyer et al. 1998; Serjeant et al. 2002). 

Optical spectroscopy is available in Gallego et al. (1996), and photometry in the Gunn-r 
and Johnson-5 filters can be found in Vitores et al. (1996) and Perez-Gonzalez et al. (2000), 
respectively. NIR data (J and K s bands) were presented in Gil de Paz et al. (2000) and Perez- 
Gonzalez et al. (2003a). An Ha imaging study has also been carried out recently (Perez- 
Gonzalez et al. 2003). For a complete compendium of the spectroscopic and photometric 
properties of the sample, the reader may refer to Perez- Gonzalez et al. (2003a, and references 
therein) and http:/ /www. ucm.es/info/Astrof/UCM_Survey/survey.html. 

In this Letter we present LFs in the 4 available optical and NIR bands: Johnson-5, 
Gunn-r, J and K s . Out of the original 191 galaxies, 15 were classified as Active Galactic 
Nuclei by Gallego et al. (1996) and have been excluded. A further 11 galaxies have no data 
for the line+continuum magnitudes (see Section 3). Thus, the sample used in this Letter is 
composed of 165 objects. The total numbers of objects with available photometry in each 
band are: 165 in B, 142 in r, 143 in J and 164 in K s . 

The luminosities have been corrected for Galactic extinction using the maps of Schlegel 
et al. (1998) and the extinction curve in Cardelli et al. (1989). We have also applied k- 
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corrections from Fioc & Rocca-Volmerange (1999) for BJK S and Fukugita et al. (1995) for 
Gunn r, taking into account the morphological types of the UCM galaxies (Perez-Gonzalez 
et al. 2001). The k-corrections applied are small because of the low redshifts of our galaxies 
(z < 0.045). 

Total stellar masses for nearly all the UCM galaxies (156 objects) were calculated in 
Perez-Gonzalez et al. (2003b) using a complete set of stellar population synthesis models to 
reproduce the broad-band and emission-line spectral features in the optical and NIR. Taking 
into account different star formation histories, the technique calculated mass-to-light ratios 
in the i^-band (A4/Lk s ) on a galaxy-by-galaxy basis. This Letter uses the results on masses 
and burst strengths (ratio of the mass of the stars younger than 10 Myr and the total stellar 
mass) which were calculated in Perez-Gonzalez et al. (2003b), assuming an instantaneous 
burst of star formation with a Salpeter (1955) initial mass function (M\ ow = 0.1 .M© and 
M up = 100 .M©) occurring on a 'normal' spiral galaxy. We used the stellar evolution syn- 
thesis library of Bruzual & Chariot (1999, private communication). The attenuation of the 
burst luminosity was modelled with the recipe given in Chariot & Fall (2000). The stellar 
masses obtained using the alternative extinction recipes considered in Perez-Gonzalez et al. 
(2003b) were statistically very similar. Thus, the choice of extinction recipe has negligible 
effect on the results presented here. 



3. Estimating the luminosity and stellar mass functions 

The selection of the UCM Survey galaxies depends primarily on the Ha line flux and its 
contrast against the continuum (Gallego et al. 1995). This contrast may be measured by a 
line+continuum magnitude m^c (Salzer 1989). In the present work we are concerned with 
the calculation of the luminosity and mass functions of the sample in broad-band filters far 
away from the passband were the selection took place. Therefore, a proper way of estimating 
these functions is using a bivariate luminosity function (BLF), <fi(M L+c , M n ), which allows 
for sample selection in one particular magnitude (m L+c ), and may be integrated over this 
magnitude to obtain the LF in the required filter, <f>(M n ). The method is analogous to the 
one described in Loveday (2000). It is an extension of the stepwise maximum-likelihood 
technique presented in Efstathiou et al. (1988), which properly handles the effects of not 
observing a small fraction of the entire sample in the band where we are estimating the LF 
(Sodre & Lahav 1993). 

All the LF data points have been fitted to Schechter (1976) functions taking the obser- 
vational errors into account. The uncertainties in the fitted parameters have been estimated 
considering uncorrelated random variations of the LF points following a Gaussian distribu- 
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tion with its width set by the errors of the data points. We calculate the errors in a, M* 
and 0* by repeating the fit for each set of points with small random variations. 

The bivariate method of estimating LFs has also been applied to the calculation of the 
total stellar mass and burst mass functions of the UCM sample, i.e., the volume density of 
local star-forming galaxies and local starbursts as a function of mass. The use of this proce- 
dure when estimating the SMF relies on principles similar to those of the LF calculations. It 
uses a bivariate function <J)(Ml + c, Ai*) defined as the number density of galaxies of a given 
Ml + c and M.*. Note that the SMF is often obtained by multiplying the i^ s -band LF by 
a constant M./Lk s - However, several authors (Brinchmann & Ellis 2000; Bell & de Jong 
2001; Papovich et al. 2001; Perez-Gonzalez et al. 2003b) have suggested there could be size- 
able galaxy-to-galaxy variations in M./ L Ks because of differences in star formation histories. 
Therefore, it seems safer to estimate the SMF using individually-estimated M/L Ka ratios 
for each galaxy, as we do here. 



4. Results 

4.1. Luminosity functions 

In order to test the BLF method, we used it to re-calculate the Ha LF of the UCM 
sample. The results are almost identical to the ones achieved with the V/V max method 
(Schmidt 1968; Huchra & Sargent 1973) used by Gallego et al. (1995) and Perez- Gonzalez 
et al. (2003). This adds confidence to our results. 

Figure 1 shows the LFs of the UCM sample in BrJK s . For comparison, we also plot the 
LFs of the local galaxy population obtained from different surveys. The dotted line in the 
left panel of Figure 1 shows the S-band LF of the 2dF Galaxy Redshift Survey (2dFGRS, 
Norberg et al. 2002). These galaxies have redshifts up to ~ 0.2, with a modal z of 0.05. 
The main different between the 2dFGRS LF and the UCM LF is a lower number density 
normalization (^ucm/^mfgrs = 0-19)- This is not a surprise since not all the galaxies in 
the 2dFGRS emit in Ha. The r-band LF (second panel of Figure 1) is compared to the 
LFs derived from the Las Campanas Redshift Survey (LCRS; dash-dot line; (z) ~ 0.1; Lin 
et al. 1996), and the Century Survey (CS; dashed line; z < 0.15; Geller et al. 1997) 5 . The 
comparison reveals again a lower number density of UCM galaxies (0ucm/0lcrs = 0-15 



5 These LFs have been transformed into the r-band using r—Rc = 0.36, typical of spiral galaxies (Fukugita 
et al. 1995). Note also that these LFs have been obtained using different Q,m and A values, but given the 
moderate redshift of the samples we have not made any correction in order to match the cosmologies. 



- 6- 



and 0ucm/0cs = 0.11). The J- and fC-band LFs are compared with those of the 2dFGRS 
(dotted lines in the last two panels of Figure 1; Cole et al. 2001), and the Two Micron All 
Sky Survey (2MASS; dashed lines; Jarrett et al. 2000). The J-band 2MASS LF comes from 
Balogh et al. (2001), with the global normalization of Cole et al. (2001). The K s 2MASS LF 
was published by Kochanek et al. (2001). In both bands a normalization offset is observed 
again: 0u CM /02dFGRS = °- 27 in J ; 0ucmMc1fgrs = °- 26 and 0ucmMmass = °- 24 in K s- 

The shapes of the luminosity distributions of the UCM sample are not very different 
from those of the global population of galaxies in the local Universe. However, since the 
UCM galaxies have typically higher star formation activity than "normal" quiescent spiral 
galaxies (Perez-Gonzalez et al. 2003b), a lower number density is found. The ratios between 
the integrated luminosity densities for the local star-forming objects in the UCM Survey and 
for the global population of galaxies are 0.15 ± 0.05, 0.11 ± 0.04, 0.23 ± 0.04, 0.18 ± 0.03 
in B, r, J, and K s . Note that the local galaxy population LFs in the B, J and K s bands 
have been derived in a homogeneous way from the 2dFGRS. Remarkably, these luminosity 
density ratios are comparable at all wavelengths and very similar to the ratio of the number of 
enhanced star- forming galaxies to the total (Zamorano et al. 1996). The luminosity density 
ratios do not depend strongly on luminosity. Moreover, these ratios are a factor of ~ 2 
higher than the nearby galaxy-pair fraction, estimated to be a 6-10% by Keel & van Soest 
(1992). This might suggest that a significant fraction of the star-formation activity in the 
local Universe is driven by minor mergers, where the smaller in-falling object is not always 
detectable. 

In the NIR, the UCM galaxies yield fainter M* values (by 0.4 mag) than the 'normal' 
galaxy samples. This is directly related to the lower typical stellar masses found for the local 
star- forming objects (Perez-Gonzalez et al. 2003b): the J and K s luminosities are strongly 
correlated with the total stellar mass, while being less influenced by recent star formation. 
However, the enhanced star formation is responsible for the UCM Survey presenting M* 
values similar to those of normal galaxies at optical wavelengths despite having smaller 
stellar masses. 



4.2. Mass functions 

Figure 2 shows the total stellar mass function (filled circles) and burst mass function 
(stars) for the UCM galaxies calculated as described in Section 3. Schechter function fits 
(solid and dotted lines respectively) are also shown. In addition, we have calculated the 
SMF multiplying the if s -band LF by a constant M./ Lk s = 0.78 (in solar units), the average 
mass-to-light ratio used in Perez- Gonzalez et al. (2003b) for the evolved stellar populations 
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within the UCM galaxies (open circles). At high masses, this simple estimate of the SMF 
coincides with the one obtained using individual M./L Ks values. However, large differences 
appear for low mass galaxies: these objects harbor the most intense bursts (Perez-Gonzalez 
et al. 2003b) which have a major influence in the Ai/Lx s ratios. The errors introduced by 
using a constant M. / Lk s are expected to be much larger at high-z where the burst strengths 
are higher. 

A Schechter function seems to provide a reasonable fit to the total SMF. The typical 
stellar mass {M.*) of the local star-forming galaxies in the UCM sample is a factor ~ 2 
smaller than that of the local galaxy population (Cole et al. 2001). This is directly related 
to the fainter NIR M* found for the UCM galaxies. The faint end slope (a) of our SMF and 
the one published in Cole et al. (2001) are virtually identical. 

We have integrated the total SMF in Figure 2 to obtain the total stellar mass density 
in local star-forming galaxies, p^ LG (z = 0) = io 7 - 83±0 - 07 J\4@ Mpc -3 . This corresponds to 
(13 ± 3)% of the baryon mass density in the form of stars and their remnants estimated by 
Persic & Salucci (1992), Fukugita et al. (1998) and Cole et al. (2001) for the local Universe. 

Integrating the burst SMF allows us to calculate the volume-averaged burst strength b of 
these star-forming galaxies, defined as the ratio of the burst mass density to the total stellar 
mass density. We find b(z = 0) = 0.04±0.01. This means that ~ 4% of the total stellar mass 
in local star-forming galaxies corresponds to stars younger than ~ 10 Myr (Perez-Gonzalez 
et al. 2003b). Note that the galaxies within the UCM Survey have EW(Ra) > 20 A (Gallego 
et al. 1995), so this value must be considered as a lower limit to the relative fraction of young 
stars in our local Universe. 

It is interesting to compare our measure of pf LG (z = 0) with high redshift estimates. 
Our value is 1.6 times lower than the co moving stellar mass density in spiral galaxies at 
0.4 < z < 1 (Brinchmann & Ellis 2000). Also, when compared with the z > galaxy 
samples presented in Dickinson et al. (2002), our results point towards a decline in pf LG 
from z ~ 1 to z = 0. This is similar to the behavior observed for the SFR density. One 
should be cautious when making comparisons with z > 1, since the high-z samples are 
dominated by intense star-formation and may present selection biases at least as important 
as the ones in our sample, leading to an underestimate of the global p*. 

5. Summary and conclusions 

In this Letter we have presented the luminosity and mass functions of the local star- 
forming galaxies in the UCM Survey. Given that the UCM galaxies were selected by their 
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Fig. 1. — Luminosity functions for the UCM sample in the Johnson- £?, Gunn-r, J and K s 
bands. The Schechter fit parameters are given inside each panel. The units of <p* are Mpc~ 3 . 
The LFs of several samples of 'normal' galaxies have also been plotted for comparison (see 
text for details). 




lo g(^stars) 

Fig. 2. — Total stellar mass function (filled circles and solid line) and burst mass function 
(stars and dotted line) for the UCM sample. Open symbols show the mass function calculated 
by multiplying the i^-band LF by a constant M/L Ka = 0.78. The Schechter fit parameters 
are given for the total masses (upper-left corner) and burst masses (lower-right corner). 
Masses are in solar units and 0* in Mpc~ 3 . 
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active star formation, our results are directly comparable to those of high redshift star- 
forming galaxies. In estimating these functions we have accounted for the Ha selection of 
our survey and calculated stellar masses on a galaxy-by-galaxy basis. 

The main results of our study are: (1) The luminosity functions of the UCM Survey 
galaxies in the optical and near-infrared have lower normalizations than the LFs of the 
global population of galaxies, indicating a lower number density. This is not surprising since 
the UCM galaxies have significantly higher star formation activity than 'normal' quiescent 
galaxies, and galaxies without star formation are not present in the sample. (2) In the NIR, 
UCM galaxies have fainter characteristic magnitudes (M*) than 'normal' galaxies. This 
is a direct consequence of the lower typical stellar masses of our objects. However, the 
optical M* values of the UCM galaxies are similar to those of 'normal' galaxies, suggesting 
the enhanced star formation has boosted up their optical luminosities. (3) The fraction 
of the total luminosity density in the local Universe associated with star-forming galaxies 
is 0.15 ± 0.05, 0.11 ± 0.04, 0.23 ± 0.04, 0.18 ± 0.03 in B, r, J, and K s . Remarkably, 
these luminosity fractions are comparable at all wavelengths and do not depend strongly 
on luminosity. This may be due to selection effects shared by all magnitude (and surface 
brightness) limited surveys. (4) If we fit the total stellar mass function using a Schechter 
parameterization, we obtain a = —1.15 ± 0.15, log(.M*) = 10.82 ± 0.17 M® and log(0*) = 
—3.04 ± 0.20 Mpc -3 . This gives an integrated total stellar mass density in star-forming 
galaxies in the local Universe of io 7 - 83±0 - 07 M.® Mpc -3 . (5) The volume-averaged burst 
strength of the UCM galaxies is b = 0.04 ± 0.01, defined as the ratio of the mass density 
of stars formed in recent bursts (age < 10 Myr) to the total stellar mass density in UCM 
galaxies. (6) We estimate that (13 ± 3)% of the total baryon mass density in the form of 
stars is associated with star-forming galaxies in the local Universe. 
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